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A method to generate new classes of random matrix ensembles is proposed. Random matrices 
from these ensembles are Lax matrices of classically integrable systems with a certain distribution 
of momenta and coordinates. The existence of an integrable structure permits to calculate the joint 
distribution of eigenvalues for these matrices analytically. Spectral statistics of these ensembles are 
quite unusual and in many cases give rigorously new examples of intermediate statistics. 
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Introduction: Statistical properties of surprisingly 
many different problems can be described by two main 
distributions: the Poisson statistics of independent ran- 
dom variables and the Random Matrix Theory (RMT) 
statistics. Classically integrable systems display spectral 
statistics which are close to Poisson distribution [l] but 
classically chaotic systems, in general, are well described 
by RMT 0. The same classes of spectral statistics ap- 
pear in the investigation of disordered systems, in partic- 
ular in the study of the 3-dimensional Anderson model 
(see e.g. Q and references therein). Though the univer- 
sal character of the Poisson and RMT statistics is well 
established (but not always well understood) these two 
distributions do not exhaust all possible behaviors, even 
in the 3-d Anderson model. A new phenomenon appears 
in this model when the disorder strength is set at a special 
value corresponding to the metal-insulator transition. In 
Q it was established numerically that in this case spec- 
tral statistics has features intermediate between Poisson 
and RMT behaviors. Similar hybrid statistics has been 
observed numerically Q for certain dynamical systems 
which are neither integrable nor chaotic but belong to 
the class of pseudo- integrable systems [1] . 

These and other examples demonstrate the possible 
existence of a not yet well-defined class of intermedi- 
ate spectral statistics Q with two characteristic features: 
level repulsion at small spacings as in RMT and exponen- 
tial decrease of the nearest-neighbor distribution at large 
spacings as in the Poisson distribution. The usual ran- 
dom matrix ensembles are chosen in such a way that their 
measure is invariant under conjugation M — > U Ai 
over a group of either unitary, orthogonal, or symplcctic 
matrices Q ■ The invariance of eigenvalues of M under 
these groups permits to find the exact joint distribution 
of all eigenvalues Xj. In the simplest setting 01 the dis- 
tribution reads 



P(A) - exp 



-a^A2.+/3^1n|A,-A, 



j<k 



where (3 = 1,2,4 for, respectively, orthogonal, unitary 
and symplectic ensembles. 

For intermediate statistics the situation is different. 
Physical problems giving rise to intermediate statistics 
have a natural basis in which they are defined, and in 
general do not possess explicit invariant measure, which 
makes the progress of their analytical treatment difficult. 

In this letter we introduce new families of random 
matrix ensembles which are not invariant over geomet- 
rical transformations, but still allow to obtain an exact 
joint distribution of eigenvalues analogous to H]). These 
ensembles give new non-trivial examples of intermediate 
statistics. 

General construction: To define our random matrix en- 
sembles we consider a classical one-dimensional A'^-body 
integrable model such that the equations of motion are 
equivalent to the matrix equation 



L^ML- LM 



(2) 



for a pair of Lax matrices L and M depending on mo- 
menta p and coordinates q @ . 

It is the Lax matrix L = L(p,q) that we propose to 
consider as a random matrix depending on random vari- 
ables p and q distributed according to a certain " natural" 
measure 



dL = P{p , q ) dp dq 



(3) 



which depends on the system. The only information we 
shall use from the integrability of the underlying classical 
system is the existence and explicit form of action-angle 
variables Ia{p,(l) and (/)Q,(p,q), and the identity 



n '^pj = n 



(4) 



(1) 



due to the canonicity of the action-angle transformation. 
Direct proof of this key identity is difficult and implicit 
methods were used to establish it [lo|, |Tl| . Action vari- 
ables turn out to be usually the eigenvalues Aq, of the 



2 



Lax matrix or a simple function of them. The canonical 
change of variables in ([3]) from momenta and coordinates 
to action-angle variables leads to a formal relation 



dL = 7'(A,0)dAd0 



(5) 



The exact joint distribution of eigenvalues is then ob- 
tained by integration over angle variables, which can eas- 
ily be performed in all cases considered: 



P(A) = Jv{X,cl) 



)d0 



(6) 



This scheme is general and can be adapted to several 
different models. Due to space restrictions we consider 
here only two representative ensembles, based on the 
rational Calogero-Moser (CM) [13and the trigono- 
metric Ruijsenaars-Schneider (RS) 'l3'| models. Other 
examples and details of the calculations will be presented 
elsewhere 14| . 



Calogero-Moser ensemble: The Hamiltonian of the ra- 
tional CM model reads 



k Kj 

with the following N x N Hermitian Lax matrix 
Ljk = PjOjk + i.g- 



(7) 



(8) 



Let Aq and Uk (a) be eigenvalues and right eigenvectors of 
L. In [10] it is proved that the matrix Q (called conjugate 
to L), defined by 



k 



can be written as 



1 - Sais 



(9) 



(10) 



where the new variables (jja — Qaa are angle variables 
canonically conjugated to the action variables Aq. Equa- 
tion (|10p is similar to ([5]) with the substitution g — > —g, 
Qj — !■ eigenvalues Aq. and pj — > angle variables (j)a ■ 

Consider now an ensemble of random matrices ([5]) de- 
fined by random variables pi and qj with measure 



dL exp 



-aTrL'-bJ2<ll 

k 



dp dq 



(11) 



where a and b are positive constants. Using the described 
canonical change of variables and taking into account 
that X^fc'Zfe = TrQ^ with Q given by (O, the measure 
in (jlip can be rewritten as 



dL - exp [-aTr - bTi Q^] dXdcj) 



exp 



E 



bg' 



(Aa — X/^y 



(12) 
dXdcj) . 



Integration over <p gives a constant and we are left with 
the following exact joint distribution of eigenvalues of the 
Lax matrix L with measure (|lip 



P(A) ~ exp 



1 



(A„ - Xfjy 



(13) 



According to this formula, eigenvalues of the above CM 
ensemble behave as a 1-d gas of particles with inverse 
square inter-particle potential. No long-range interaction 
proportional to InjAi — Aj| is present, in contrast with 
standard random matrix ensembles ([T]). 

The fast decrease of inter-particle potential in (fT3l) 
with the distance between particles permits to approx- 
imate (see e.g. Q) the nearest-neighbor distribution of 
eigenvalues of the Lax matrix ^ by the formula 

P{s) « Ac-^^''-^' (14) 

where B is a fitting constant, and constants A and C are 
determined from the normalization conditions. This ex- 
pression is not exact but may be considered as an analog 
of the Wigner surmise in RMT [?] . 

The measure dTTI) for coordinates q corresponds to N 
particles with repulsion confined in an interval of the or- 
der of 1/Vb. In order to simplify numerical investigation 
it is therefore natural to use the "picket fence" approxi- 
mation Qk ^ k with integer k. We thus replace the matrix 
L by a simpler matrix 

1 - Sjk 



Ljk = PkS 



•]k + ig- 



(15) 



'2(i-fc)- 

The decrease as \j — k\^^ of non-diagonal elements in 
L is a characteristic feature of intermediate systems 0]. 
In numerical calculations we chose N variables pk as 
i.i.d. random variables with uniform distribution be- 
tween — 1 and 1. In Fig. [T] we show the nearest-neighbor 
spacing distribution for the matrix L for several values of 
the parameter g. The simple surmise (jl4p is practically 
indistinguishable from numerical results (see inset), 
which confirms the existence of unusual exponentially 
strong level repulsion in this model. 

Ruijsenaars-Schneider ensemble: The second example 
we consider here is the trigonometric RS model [l^. It 
appears that the random matrix ensemble that was pro- 
posed in [l3| as a quantization of a pseudo-integrable 



17l| is a 



interval-exchange map and investigated in [H 
particular case of this model. 

The RS model is determined by the Hamiltonian 

N 

i?(p,q)=Ecos(p,)V^//'(T,q) 
where Vj (t, q) depends on q and a real parameter r as 



(16) 



^.(^,q) = n (i- 



sin^ T 



sin [{qj - qk)/2] 



(17) 
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FIG. 1: Nearest-neighbor distribution for the CM random 
matrices psp with g — 0.1, 0.5, 1, and 2 (from left to right), 
averaged over 5000 reahzations of matrices of size A'' = 301. 
Black solid lines are numerical results, red dashed hnes indi- 
cate the fit (|14p with fitted values of B. For the considered 
values oi g, B = 0.096, 0.618, 1.46, and 3.11 respectively. In- 
set: Difference between numerical result and fit for g = Q.l 
(black), 0.5 (red), 1 (green), and 2 (blue). 



The Lax matrix i(p, q) for this model is a unitary matrix 
given by 



L,-fe(p,q) = e'^(^-i)+'P^Qfe(q)e'('--«^- 



)/2 



(18) 



(we choose a phase factor different from the one in [ll| 
in order to get Eq. ^ when qu = 2nk/N). Here C(q) 
is the orthogonal matrix 



^ , X „rl/2, ^ SUIT 

sm[{qj -qk)/2 + ' 
with Wj (r, q)Wj (— T, q) = V-j (r, q) and 



1/2/ 



'T, q) 



(19) 



w^.(-,q)-n^^^S^-^^^- (20) 



sin[(qj - qk)/2] 



Action-angle variables are obtained similarly as for the 
CM model. Here one considers llj] the conjugate matrix 
Q defined by 



Qa/3=^<(«)e'«'=Ufc(/3) 



(21) 



where Uk (a) are eigenvectors of the Lax matrix (|18p cor- 
responding to eigenvalues Aq = e'^°. In ^ it is shown 
that Qa/i can be written in the form ^5)) with the fol- 
lowing substitutions: t — )■ — t, q„i action variables 9a 
and pk — > angle variables (f>a canonically conjugated to 

9a. 

The important difference of this model from e.g. the 
above CM model is that the Hamiltonian (fTB|) and the 
Lax matrix are defined not on the whole q-space 



but only on a subset of it where all Vj{a,q) in ([17]) are 
positive (notice the square roots in these expressions). 
These restrictions depend only on coordinates and on r 
(in [m only the case < r < tt/N had been considered). 
Let R{t, q) be the characteristic function of this subset 



R{T,q) 



1 when V, (r, q) > , j = 1, , 



.N 







otherwise 



(22) 



We choose as a "natural" measure for the RS ensemble 
the uniform measure of random variables p and q on the 
region allowed by the above restrictions. This implies 
that the measure on the RS ensemble is chosen as 



dL ~ R{t, q) dp dq 



(23) 



We transform this expression to action-angle variables, 
and perform the integration over angle variables. Since 
Q and L have the same form but with A and q inter- 
changed, A is subject to the same restrictions as q. We 
conclude that the exact joint probability of eigenvalues 
of the ensemble of random RS matrices (fTSl) is 



(24) 



As mentioned 
vestigated in 



is a generalization of the model in- 



and [17 1 



The simplest non-trivial 
new case corresponds to the choice r = irb/N with 
fixed b. To find R{t, A) in this case we notice that as 
in (16j and (17j the matrix (jlSp permits two rank-one 
deformations with known eigenvectors and eigenvalues 
N^j t^ = Ljte^'^^^"'?''"''^'^^. Generalizing the discussion in 
[l6| and [l3|, one can prove [l3| that for N large enough, 
there exist exactly n = [b] other eigenvalues at angular 
distance 2'Kh/N from any eigenvalue 9a (here \b] denotes 
the integer part of h) . 

Consider an ordered sequence of eigenphases on the 
unit circle, Oi < O2 < ... < 9^ and denote the nearest 
differences by — 9k+i ~ 9k- Introducing two functions 



1 when < X < 6 
otherwise 



'i{x) = l- f{x), (25) 



one can show that these restrictions give rise to the 
following expression for the joint probability ([24]) of RS 
Lax matrix eigenphases inside an interval of length A 



N 



N 



(26) 



k=l 



where sj = £,j + . . . + ^j+n-i and n = [b]. This formula 



means that eigenvalues of RS random matrices ([T8|) with 
r = Trb/N behave exactly as a 1-d gas where each particle 
interacts with n = [b] nearest-neighbors. In is known that 
in this case all correlation functions in the limit of large 
N can be calculated by the transfer operator method 
(see e.g. Q). Here we present a few results for the k^^ 
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FIG. 2: Nearest-neighbor distributions (from left to right) 
P(l,s) (black), P(2,s) (red), and P(3, s) (blue) for RS ran- 
dom matrices (|29|l with r = ivb/N and & = 9/4. All curves 
correspond to A'' = 701 averaged over 1000 realizations of ran- 
dom phases p^. Solid lines are theoretical predictions. Dashed 
vertical line indicates abscissa 9/4. 



nearest-neighbor distributions, P(k, s) which determine 
the probability that in the interval of length s there exist 
exactly fc — 1 other eigenvalues (P(l, s) = P{s)). The 
details will be discussed elsewhere [14| . 

When < 6 < 1, P{k, s) = ioi < s < kb and for 
s > kb P{k, s) is a shifted Poisson distribution 



g-(s-fcb)/(i-b) 



{s-khf 



(27) 



For larger 6 formulas, though explicit, become tedious. 
For example, for b = 4/3, P(l,s) is non-zero only when 
< s < 4/3 and P(2,s) when 4/3 < s < 8/3. Inside 
these intervals 



P{s) 



81 
64' 



16 



512 



(28) 

To simplify numerical calculations we use (as in (jlSp ) the 
picket fence approximation of coordinates qu = li^kjl^ . 
At these values of q, Wj(r, q) = sinA^r/(A^sinr) and 
the Lax matrix takes the form 



1 



^2iTAf 



7V(1 - e2'ri(j-'=)/^'+2ir) 



(29) 



with random phases "Pk uniformly distributed in [0, 27r[. 

When T = 7ra with fixed a, this matrix up to notations 
coincides with the one proposed in [l5!| and investigated 
in [lij and 17 1. We put r = ■nb/N and perform numer- 



ical calculations of spectral statistics of matrix ([29)1 for 
different values of h and find that all above formulas very 
well agree with numerics. As an illustration we present 
at Fig. [2] a case with h = 9/4 for which explicit formulas 



are too long to be presented here. Even in this more com- 
plicated case analytical results are difficult to distinguish 
from numerics. 

Conclusion: To summarize, we proposed a general 
method of constructing non-invariant random matrix en- 
sembles whose joint distribution of eigenvalues can be 
calculated analytically. These ensembles are Lax ma- 
trices of classically integrable TV-body models, equipped 
with a suitably chosen measure of momenta and coor- 
dinates which depends on the model. For such matrix 
ensembles the symmetry groups of usual RMT are re- 
placed by the underlying structure of integrable flows 
generated by N conserved quantities. It is this struc- 
ture which makes possible the explicit construction of 
joint probability of eigenvalues in these ensembles. Spec- 
tral statistics of these ensembles are quite unusual and in 
many cases they present new examples of non-universal 
intermediate statistics. In all considered cases eigenfunc- 
tions computed numerically present multifractal proper- 
ties , which is a typical feature of intermediate statis- 
tics [3|. It is interesting to note that a specific random 
matrix ensemble, which appeared in [l^ as the result of 
quantization of an interval-exchange map, belongs to this 
class. 
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